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ABSTRACT 



A computer simulation of an aircraft making a low level penetra- 
tion over hostile terrain is described. The model discussed is intended 
as a tool for the determination of appropriate defensive tactics for air- 
craft upon receiving fire from conventional ground weapons. In any com- 
puter model of a ground to air combat environment it is essential to have 
a realistic simulation of the terrain. The method used in this simulation 
is discussed and shown to be extremely accurate. 
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1 . 



Introduction 



PENAIR is a time step computer simulation of an aircraft flight over 
hostile terrain. It is an analytic computer war game in which one aircraft 
and a maximum of two hundred weapons are the principles . Any piece of 
actual terrain may be simulated using a contour map of the area to obtain 
the input parameters . 

The simulation is intended as a tool for the determination of appropri- 
ate defensive tactics for an aircraft upon receiving fire from conventional 
ground weapons. Using the characteristics of the aircraft and of the wea- 
pons the model determines the survival probability of the aircraft on any 
of several flightpaths . There are no restrictions on the aircraft's speed, 
altitude or other characteristics and therefore this model is equally suita- 
ble for helicopters or jets. 

The aircraft flightpath can be predetermined or it can be generated 
by the model. Upon receiving fire the aircraft can be programmed to exe- 
cute any of 18 evasive maneuvers . A comparison of the survival probabili- 
ties can then be made thus giving an indication of the appropriate defensive 
maneuver to employ. 

The terrain simulation and generated flightpaths appear to be ex- 
tremely realistic. The terrain is simulated by fitting, in the least squares 
sense, a polynomial function of two independent variables (x,y) to values 
of a dependent variable specified at points on a rectangular grid in the 
plane of the independent variables . This is accomplished by a separate 
computer program, TERRAIN, which generates the necessary inputs for 
PENAIR. 

2 . Design of the Model 

PENAIR is a computer simulation of a single aircraft making a low 
level penetration over terrain containing hostile small arms and con- 
ventional anti-aircraft weapons. Any actual or fictitious terrain may be 
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simulated. No offensive aircraft weapons are considered although 18 dif- 
ferent defensive maneuvers such as increasing speed and diving for the 
deck may be performed. Ten types of ground weapons may be used with 
a maximum of 20 weapons per type. The program computes the aircraft's 
time in view of each weapon, the probability of the aircraft surviving each 
weapon and the aircraft's total survival probability. The simulation runs 
in discrete time steps using any time interval desired. The program does 
not utilize Monte Carlo methods; aircraft survival probability for a speci- 
fied flightpath over a specified terrain is computed in a single run . Flow 
chart 1.0 of Appendix I f page 30 , indicates the general logic of the PENAIR 
program . 

PENAIR is written in FORTRAN-63 for the CDC 1604 computer. The 
PENAIR program is included as Appendix II. A separate program, TERRAIN, 
is utilized to supply inputs to PENAIR for terrain simulation. The terrain 
input program supplies the necessary data on punched cards for the PENAIR 
program. The TERRAIN program is also programmed in FORTRAN-63 for the 
CDC 1604 and is included as Appendix IV. 

The terrain used in the PENAIR program may be divided into as many 
as 12 separate squares for simulation. Each of the 12 individual squares 
is approximated by a separate polynomial. The 12 squares may be of any 
common size. Figure 1 illustrates the arrangement of these squares. The 
coordinate system has been established such that the lower left comer is 
the origin with increasing x to the right and increasing y upward. It is 
possible to use actual military grid coordinates since the origin may 
take any desired x and y values . 

A special PENAIR output is a set of punched cards to be used as 
input to a graph plot routine for an off-line plotter. This routine then draws 
the actual flightpath over the simulated terrain. Output from this routine 
is included in Section 6. 
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THE ARRANGEMENT OF THE INDIVIDUAL 
SQUARES OF TERRAIN SIMULATED IN 
THE PENAIR PROGRAM 



3 . Aircraft Movement 

The PENAIR model has provisions for three types of simulated aircraft 
flightpaths. The first of these flightpaths allow the PENAIR program to 
select one of 18 different evasive maneuvers for the aircraft upon receiv- 
ing fire. The 18 evasive maneuvers consist of six basic maneuvers with 
three possible speed changes for each. The six basic maneuvers are: 

(1) continue flying straight ahead at the same altitude, (2) commence a 
climb straight ahead at the aircraft's maximum rate of ascent, (3) turn 
away from the weapon remaining at the same altitude, (4) turn away and 
commence a maximum rate climb, (5) dive for the deck and commence nap 
of the earth (contour) flight straight ahead, or (6) turn away from the 
weapon, dive for the deck and commence nap of the earth flight. For each 
of the above maneuvers the aircraft may decrease speed, remain at the 
same speed or increase speed. All of the above maneuvers are gener- 
ated using the input characteristics of the aircraft being simulated. It 
should be realized by the PENAIR user however that only one basic maneuver 
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may be employed in one computer run. The maneuver employed is deter- 
mined by the numerical value of the input flag IMARK. 

The other two simulated aircraft flightpaths are preplanned. In the 
first of these, the flightpath is completely determined in advance and 
read into the PENAIR program from data cards. In the second, the x and y 
coordinates of points along the flightpath are read in to the PENAIR pro- 
gram. From these points the PENAIR program computes a nap of the earth 
flightpath for the aircraft over the simulated terrain. In both of these 
cases, the PENAIR program makes no deviation from the pre-planned flight- 
path and the aircraft makes no use of the evasive maneuvers. A complete 
description of the inputs needed to regulate aircraft movement is included 
in Appendix III. 

4. Determination of the Aircraft’s Survival Probability. 

The major output of the PENAIR program is the probability of aircraft 
survival . This survival probability is computed as a function of weapon 
type, range and time in view. The necessary inputs are weapon character- 
istics and the single time-interval kill probabilities for each weapon. 

Given these kill probabilities, the probability of an aircraft surviving any 
single weapon is computed as: 

P (Survival) = (1 - P ) (1- P ) . . .(1 -p ) 

i £ n 

where P^ is the probability of kill for one time interval at range i for a 
given weapon. The total probability of aircraft survival is then computed 
as the product of the individual probabilities of survival for each weapon. 
(See Flow Chart 5.1, Appendix I, page 48 ) . It is assumed that the proba- 
bility of a kill by one weapon is independent of all other weapons . The 
number of time intervals the aircraft receives fire depends on the time in 
view, minimum and maximum weapon range, and the weapon acquisition, 
firing , and reloading times . 
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Time in view is the one parameter to be determined. The method used 
by this simulation to determine time in view will be explained using Figure 
2 and Flow Chart 7.0 of Appendix I, page 49. 

The aircraft is considered in view of a weapon if every point on the 
line of sight WA between the weapon position W, and the aircraft A, is 
higher than the corresponding point of the terrain. To determine this the 
x,y coordinates of the point on the line TR, 10 meters from the weapon are 
used to calculate the elevation of the terrain at that point. This elevation 
is then compared with the height of the line of sight at those coordinates . 

If the terrain is lower, another point 50 meters farther along the line TR 
is chosen and the same comparison is made. If the terrain elevation is 
below the line of sight elevation at this point this procedure is repeated 
every 50 meters. This process is continued until it is determined that the 
aircraft is in view of the weapon or that the terrain elevation exceeds the 
line of sight. If the latter occurs the aircraft is considered not in view 
of this weapon at this time. This process is repeated for all weapons and 
then the aircraft is advanced one time interval. The number of consecutive 
time intervals the aircraft remains in view of any one weapon is then used 
to compute the current total amount of time in view per weapon. 




DETERMINATION OF THE AIRCRAFT'S TIME 
IN VIEW OF A WEAPON 
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Two other methods of determining if the aircraft was in view of the 
weapon were tested. Elevation comparisons between line of sight and 
terrain were made at 25 meter increments on the line of sight WA, and at 
25 meter increments along TR. No significant differences were noted in 
the times of view computed by the 3 methods mentioned. However, the 
method outlined in the previous paragraph greatly reduces computer run- 
ning time for the PENAIR program . 

After the aircraft has been in view of a weapon for the required acqui- 
sition time, the weapon minimum and maximum ranges are checked. If 
the aircraft is within firing range, the weapon commences firing and the 
probability of survival is determined for this weapon in this time interval. 
The probabilities of aircraft survival for each weapon are then computed 
for all following time intervals until the aircraft either goes out of view 
and/or range or until the weapon has expended its ammunition. In both 
cases no further probabilities are computed until the weapon is reloaded 
or until the aircraft comes back in view and within range. It is assumed 
that after reloading reacquisition requires one-half the acquisition time. 

If the aircraft goes out of view the weapon must wait its full acquisition 
time after the aircraft is again in view. 

5 . Terrain Simulation . 

In any computer model of a ground combat environment it is essential 
to have a realistic simulation of the terrain in order to determine the 
mask angles of ground weapon positions. To do this the simulation must 
provide the elevation of any point in the combat area. Present computer 
methods of terrain simulation either use very specialized equipment or 
involve prohibitive amounts of computer storage. The common method is 
to divide the terrain of interest into equal areas and assign an appropriate 
elevation to each. The assigned elevation might be the average elevation 
of the area. In this method it is assumed that all points in a given area, 
for the purpose of simulation, have the same elevation. This is a 



o 



crude approximation of terrain features for large areas. The approximation 
approaches reality as the size of the areas decreases. However, since 
the number of the areas considered is inversely proportional to the size of 
the areas, as the approximation of the terrain approaches reality the com- 
puter storage required becomes very large. 

For the model developed here the terrain features are expressed in 
polynomial form, z=f(x,y) . The general problem is that of fitting points 
in three dimensions to a surface represented by the polynomial, z = f(x,y) . 

A method of forming this polynomial using the principle of least squares, 
given the set of data points { (x,y,z) } , had been programmed in FORTRAN-60 
for the CDC 1604 [21 . This program was converted to FORTRAN-63 and a suc- 
cessful attempt was made in fitting a polynomial to a fairly rough surface 
(mountainous terrain.) The polynomials that best fit mountainous terrain 
whose elevations ranged from 800 to 1300 feet have approximately 400 terms 
with the maximum degree of both x and y near 20. 

The program designed to fit the polynomial to the terrain required the 
input points be from an equally spaced grid in the x,y plane. For this thesis , 
areas of 1000 meters on a side were arbitrarily chosen from maps of actual 
terrain. In order that the fit be acceptable on the edges of the grid square, 
an area of at least 1200 x 1200 meters with the desired 1000 meter square 
centered inside had to be considered. Consecutive x and y coordinates of 
each individual contour line were recorded and a program was written to 
produce equally spaced elevations from the consecutive data points. How- 
ever, because of the great amount of storage required, these two programs 
could only be incorporated into one by the use of overlays. An overlay is 
a feature of FORTRAN-63 designed to increase the use of the core storage 
of the computer by using magnetic tape as auxiliary storage. 

The program to compute z= f(x,y) from the consecutive (x,y) points 
of actual contour lines consists of a main program and three overlays . 
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The sole purpose of the main program is to call the individual overlays 
into operation. The first overlay, MAPPING, computes the equally 
spaced elevation points using the (x,y) points from the contour lines as 
input.. The second overlay, SURFIT, by using the points computed in the 
first overlay, computes the polynomial z= f(x,y) and records the co- 
efficients on punched cards . The second overlay also computes the ele- 
vations of a grid for use in the third overlay. The third overlay, PLOTTER, 
generates data and stores it on a magnetic tape suitable for a CDC-160 
program used to drive a CalComp 165 graph plotter. The contour map drawn 
can then easily be compared to the original map to estimate the goodness 
of fit of the polynomial . The drawing of this map is optional and is con- 
trolled by an input parameter to the main program. A detailed description 
of each overlay follows . 

Overlay 1 determines the elevation of the terrain at equally spaced 
points of a grid in the x,y plane. The program was written by T. K. Rode- 
burg of the Norwegian Defense Institute for Stanford Research Institute, 

Fort Ord, California. This overlay uses as input the lists of the x and y 
coordinates at successive points on each contour line read from an en- 
larged map of the actual terrain. The following steps are taken by this over- 
lay in the process of determining the elevation of the grid points: 

(1) The point coordinates are corrected for distortion caused by un- 
equal magnification of the original map. 

(2) Equally spaced vertical lines are determined and the coordinates 
of the intersection of the contour lines and these vertical lines are calcu- 
lated. 

(3) The elevation at these intersections are calculated by linear 
interpolation . 

(4) By further linear interpolation, elevations at equally spaced 
points on the vertical lines are found, thus obtaining the elevations of 
points on a grid in the x,y plane. 
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(5) Using the original data points, the above interpolation procedure 
is repeated at equally spaced horizontal lines. 

(6) The average of the two values thus obtained is used as the ele- 
vation at each point of the grid . 

Overlay 2 fits, in the least squares sense, a polynomial function of 
two independent variables, x and y, to values of a dependent variable, z, 
specified at points on a rectangular grid in the plane of the independent 
variables. That is, given a set of ground elevations, z at grid coordinates 

N M . m 

x and y, this overlay finds the polynomial, EE A X Y , which best 

n= o m=o nm 

fits the data. The use of orthogonal polynomials leads to a particularly 
simple system of linear equations rather than the ill-conditioned system 
which arises from the usual least squares normal equations. It also provides 
a measure of the improvements resulting from each new term included which 
further leads, in this algorithm, to an automatic selection of a "best" degree 
polynomial function as determined by Gauss' criterion [1]. It therefore de- 
termines the optimum values of N and M as well as the coefficients . 

This algorithm was programmed in FORTRAN-60 for the CDC 1604 by 
Mrs . Jean Bow of the USNPGS from an ALGOL program [2] . Minor modifi- 
cations were made to convert it to a FORTRAN-63 overlay. The FORTRAN-60 
version is a very general program whereas the overlay has had many of the 
generalities removed for ease of handling inputs which in turn reduced the 
computer running time significantly. A FORTRAN listing of this program 
appears in Appendix IV. 

The following steps are taken by Overlay 2 in the process of determin- 
ing an approximating polynomial: 

(1) the variables are normalized by subtracting off their means, 

(2) the orthogonal polynomials are evaluated , 

(3) the contribution of each orthogonal polynomial to the minimization 
of the residuals is evaluated, 



9 






,i 



(4) Gauss' criterion is applied to determine the degree of poly- 
nomial which yields the closest fit to the given data. 

(5) the orthogonal polynomial coefficients are evaluated, 

(6) the dependent variables are evaluated using the approximating 
polynomial , and 

(7) the terrain elevations utilizing a 50 meter grid are evaluated for 
input into Overlay 3 . 

Overlay 3 generates lists of x and y coordinates of successive 
points on each contour line and stores these coordinates on an output 
tape in a form suitable for a CDC-160 program used to drive a CalComp 165 
plotter. This tape is used as an input tape for the plotter which then recon- 
structs the terrain by connecting the points with straight line segments „ 

The FORTRAN-63 program used in Overlay 3 was originally developed 
for use with the IBM 7090 system and an x,y plotter, to give a complete 
traditionally styled contour map. This program was converted to FORTRAN- 
63 by Mr. R. R. Hillery of the USNPGS, for use on the CDC 1604. Subrou- 
tines used to compute the points to be plotted were provided by Dr„ M. O. 
Dayhoff of the National Biomedical Research Foundation, Silver Spring, 
Maryland [3] . Subroutines used internally to write the output tape were 
programmed by Mr. J. R. Ward of the USNPGS for use with a CalComp 
165 graph-plotting routine, DRAW. 

Figure 3 is an actual contour map from which approximately 1200 input 
points have been read. A polynomial of 400 terms was fitted and the ele- 
vations of a grid (spacing of 50 meters) were computed to plot Figure 4. 

The coefficients of this polynomial ranged from 904.18 for the constant 

-94 19 19 

term to -4.593 x 10 for the x y term. 

An analysis of a table of differences showed the maximum difference 
of 625 actual and computed elevations to be 31 feet in an elevation range 
of 900 to 1200 feet. 
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ACTUAL CONTOUR MAP OF GRID SQUARE 7174 

U. S. ARMY HUNTER-LIGGETT MILITARY 
RESERVATION 




FIGURE 4 



COMPUTER SIMULATION OF GRID SOUARE 7174 
U, So ARMY HUNTER-LIGGETT MILITARY 
RESERVATION 
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FIGURE 5 

REPRESENTATIVE FIT AT THE JUNCTION 
OF TWO GRID SQUARES 
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Figure 5 shows a representative fit at the junction of two of the grid 
squares approximated by separate polynomials. The accuracy with which 
the grid squares fit together is determined by the accuracy of the input 
data. 

6. Analysis of the Model 

PENAIR has been checked out with the use of a simple scenario. A 
3000 x 4000 meter mountainous section of the U. S. Army Military Reser- 
vation was selected for simulation. Three 10-men enemy squads armed 
with M-14 rifles were situated as illustrated in Figure 6. A UH-1B heli- 
copter was assigned the mission of delivering supplies from an airfield 
at position a to an outpost at position b „ Two test runs were made with 
altitudes of the helicopter of 1000 feet and 500 feet above the terrain. In 
both tests all 18 evasive maneuvers were checked. Some representative 
flightpaths generated by the program are illustrated in Figures 6, 7 , 8, and 
9 . The survival probabilities of the aircraft on the various attempts were 
computed and are shown in Tables 1 and 2 . All of the computed proba- 
bilities are realistic and are, in a comparative sense, what one would 
expect. 

A sensitivity analysis was then conducted on the rate of climb and des- 
cent and changes in aircraft speeds . A helicopter was simulated flying at 
100 knots with its maximum rate of climb and descent taking the values of 
500,1000 and 2000 feet per minute. The effects of these three rates on 
nap of earth flight are illustrated in Figure 10 . It should be noted that 
Figure 10 illustrates how the program generates the flightpath to overshoot 
the crest of a hill before descending just as a pilot tends to do. Figure 10 
also shows, as should be expected, the greater the rate of climb and des- 
cent, the closer the aircraft remains to the earth. Figure 11 illustrates an 
aircraft flown with a maximum rate of climb and descent of 1500 feet per 
minute and speeds of 100, 200 and 300 knots. The resulting flightpaths also 
appear realistic. 
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The flightpath computed by PENAIR illustrated in Figure 12 is of a 
helicopter making a level turn at 70, 90 and 110 knots. Here it is noted, 
for the computed flightpath, the radius of turn increases as the airspeed 
increases . 

The last analysis made was the comparison of computed results with 
actual results. Figure 13 is the simulated flightpath and terrain of an 
actual flightpath flown by a UH-1B over an actual terrain. The 44 weapon 
positions illustrated were actually manned by U „ S. Army personnel. The 
observed time the helicopter was in view of each of these positions is 
listed in Table 3 . The times in view for the simulated terrain and flight- 
path were computed and are also listed in Table 3 . It is believed that with 
better terrain data an even closer comparison could be obtained. 
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FIGURE 6 

t 

A STRAIGHT AND LEVEL FLIGHTPATH 
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FIGURE 7 

A FLIGHTPATH OF AN AIRCRAFT CLIMBING STRAIGHT AHEAD 
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FIGURE 8 

A FLIGHTPATH OF AN AIRCRAFT TURNING AWAY FROM A FIRING WEAPON 



17 




FIGURE 9 

A NAP OF THE EARTH FLIGHTPATH 
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Summary of Survival Probability of A UH-1B Helicopter Fired on by Three 
10 -Men Squads of M-14 Rifles 



Initial Speed 90 Kts Initial Altitude 650 Meters 



Evasive Maneuver 
(Commenced after receiving fire) 


Evasive 

Speed 


Survival Probabilities 
Wpn A Wpn B Wpn C 


Total 


Fly straight and level 


70 


.8576 


.8078 


1.0000 


.6928 




90 


.8723 


.8649 


1.0000 


.7545 




110 


.8723 


.8723 


1.0000 


.7610 


Climb straight ahead 


70 


.8723 


.8723 


1.0000 


.7610 




90 


.8723 


.8723 


1.0000 


.7610 




110 


.8723 


.8723 


1.0000 


.7610 


Turn away, level 


70 


.9026 


1.0000 


1.0000 


.9026 




90 


.9026 


1.0000 


1.0000 


.9026 




110 


.9026 


1.0000 


1.0000 


.9026 


Turn away and climb 


70 


.9026 


1.0000 


1.0000 


.9026 




90 


.9026 


1.0000 


1.0000 


.9026 




.110 


.9026 


1.0000 


1.0000 


.9026 


Dive and fly NOE 


70 


.8359 


.7252 


1.0000 


.6062 




90 


.8723 


.80 31 


1.0000 


.7005 




110 


.8723 


.8242 


1.0000 


.7190 


Turn away and fly NOE 


70 


.8798 


1.0000 


1.0000 


.8798 




90 


.8874 


1.0000 


1.0000 


.8874 




no 


.8874 


1.0000 


1.0000 


.8874 



Table 1 
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Summary of Survival Probabilities of UH-1B Helicopter Fired on by Three 

10-Men Squads Armed with M-14 Rifles 



Initial Speed 90 Kts Initial Altitude 450 Meters 



Evasive Maneuver 


Evasive 


Survival 


Probabilities 


(Commenced after receiving fire) 


Speed 


Wpn A 


Wpn B 


Wpn C Total 


Fly straight and level 


70 


Values not determined 






90 


at these speeds 






110 


.8309 


.7788 


1.0000 


.6471 


Climb straight ahead 


70 


Values not determined 






90 


at these speeds 






no 


.8377 


.8242 


1.0000 


.6905 


Turn away, level 


70 


.8874 


1.0000 


1.0000 


.8874 




90 


.8874 


1.0000 


1.0000 


.8874 




no 


.8874 


1.0000 


1.0000 


.8874 


Turn away and climb 


70 


.8798 


1.0000 


1.0000 


.8798 




90 


.8798 


1.0000 


1.0000 


.8798 




no 


.8798 


1.0000 


1.0000 


.8798 


Dive and fly NOE 


70 


.7558 


.9501 


1.0000 


.7181 




90 


.8027 


.9664 


1.0000 


.7758 




no 


.3306 


1.0000 


1.0000 


.8306 


Turn away and fly NOE 


70 


.8798 


1.0000 


1.0000 


.8798 




90 


.8798 


1.0000 


1.0000 


.8798 




no 


.8798 


1.0000 


1.0000 


.8798 



Table 2 




RATE OF CLIMB = RATE OF DESCENT = 500 FEET PER MINUTE 




RATE OF CLIMB = RATE OF DESCENT = 1000 FEET PER MINUTE 




RATE OF CLIMB = RATE OF DESCENT = 2000 FEET PER MINUTE 



FIGURE 10 

COMPARISON OF RATES OF CLIMBS AND DESCENTS IN A NAP 
OF THE EARTH FLIGHT AT A CONSTANT SPEED OF 100 KTS. . 
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.100 Kts 





300 KTS. 
FIGURE 11 



COMPARISON OF AIRCRAFT SPEEDS ON A NAP OF THE EARTH FLIGHT 







FIGURE 12 

COMPARISON OF THE AIRCRAFT' S RADIUS OF TURN 
AT 3 AIRSPEEDS, 70, 90 and 110 KTS. 
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FIGURE 13 



COMPUTER SIMULATION OF AN ACTUAL FLIGHTPATH 
FLOWN BY A UH=1B HELICOPTER 
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Comparison of Actual and Computed Time Intervals in View of 44 Weapon Sites 



Time Interval = .6 sec. 



Wpn. 

No 


Time 

Actual 


in View 
Computed 


Wpn 

No 


Time 

Actual 


In View 
Computed 


Wpn 

No 


Time 

Actual 


In View 
Computed 


1 


0 


0 


16 


50 


56 


31 


0 


0 


2 


0 


9 


17 


0 


0 


32 


22 


39 


3 


30 


30 


18 


43 


47 


33 


0 


0 


4 


0 


0 


19 


63 


108 


34 


89 


88 


5 


0 


0 


20 


4 


30 


35 


93 


91 


6 


0 


59 


21 


17 


74 


36 


69 


76 


7 


0 


0 


22 


0 


72 


37 


66 


58 


8 


0 


40 


23 


0 


0 


38 


49 


64 


9 


0 


0 


24 


0 


0 


39 


16 


36 


10 


0 


17 


25 


59 


63 


40 


0 


16 


11 


9 


2 


26 


0 


34 


41 


8 


8 


12 


0 


0 


27 


63 


64 


42 


56 


56 


13 


26 


39 


28 


43 


45 


43 


19 


15 


14 


40 


63 


29 


0 


9 


44 


103 
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7. Conclusions and recommendations. 

The simulation of the terrain obtained by this program has proven to be 
extremely realistic. (See Figure 4, Section 5, page 11 ) . The junctions of 
squares approximated by separate polynomials are as good as the data used 
for the approximation. (See Figure 5, Section 5, page 12 ) . This program was 
written in its present form because of the availability of input data to the 
author in the format described in Appendix V. If the actual elevations could 
be obtained using an equally spaced grid for the x and y coordinates, the 
first overlay could be eliminated It is felt that this is feasible and would 
make the preparation of the input cards much simpler. 

Among the uses of the simulation of the terrain by this method include 
ground to ground combat and ground to air combat environments. Mask angles, 
terrain profiles and actual distances between points in a mountainous terrain 
are all easily determined. The uses of the PENAIR program include the deter- 
mination of optimal flightpaths across hostile terrain, the survival proba- 
bility and optimal maneuvers to perform if the aircraft is fired upon. 

The results in Section 6 all tend to illustrate the realism of the aircraft 
penetration model. Whether the actual kill probabilities used are exact is 
not extremely important since an appropriate defensive maneuver for an un- 
armed aircraft to initiate upon receiving ground fire may be obtained by ob- 
serving the relative survival probabilities of several maneuvers. 

It is believed that with additional programming this model could be ex- 
tended to include formations of two or more aircraft, thus greatly increasing 
its usefulness. In addition, since the terrain is simulated by mathematical 
functions and since the line of sight between aircraft and weapon may be 
described by a linear relationship, it appears possible to solve the two rela- 
tions for a point of intersection. Upon accomplishing this programming, it 
is believed the computer running time would be reduced significantly. 
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APPENDIX I 

Flow Charts of the PENAIR Program 



In the flow charts of this appendix, the following conventions 
are used: 





This symbol is used for a READ 
statement. 



This symbol is used for a PRINT 
statement. 



A rectangle indicates one or more 
FORTRAN statements . A number such 
as 4.1.8 indicates the statement 
in this rectangle is further explained 
by Flow Chart 4.1.8 in this appendix. 
A number such as 420 is the statement 
number as it appears in the FORTRAN 
listing in Appendix II. 




This symbol is used for an IF statement. 




This symbol indicates the logic of this 
Flow Chart is continued at the beginning 
of Flow Chart 5.1 of this appendix. 
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10 

comma. 




This symbol indicates a connector between 
two parts of the same Flow Chart of this 
appendix. 

This symbol indicates the logic of this 
Flow Chart is continued where the same 
symbol appears in Flow Chart 4.1 of this 
appendix. 

This is the symbol for a CONTINUE state- 
ment whose statement number is 10. 

This is the symbol for a computed 
GO TO statement. 



This symbol is used for the statement 
DO 10 1= 1, 12 



BEjO/N 



B.ND 



This symbol indicates the beginning or 
entry point for the Flow Chart . 

This symbol indicates the end or exit 
point for the Flow Chart. 
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APPENDIX II 



FORTRAN Listing of the PENAIR Program 

PENAIR consists of the main program and four subroutines. The 
variable names used in the PENAIR program are defined in the comments 
appearing in the FORTRAN listing. The variable names of input and out- 
put parameters are explained in Appendix III. 
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PRINT 2041 PEN0145 

2041 FORMAT (1H1) PEN0146 

PEN0147 

COMPUTE THE MAXIMUM ANGLE OF BANK OF THE A/C. PEN014S 

AB=ACOSF( l./GMAX) PEN0149 
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PRINT 3010, TPS PEN0323 

3010 FORMAT ( ///30X ,26HTOTAL SURVIVAL PROBABILITY ,F6.A/1H1) PEN0326 

END FILE 52 PEN0327 

WRITE (1,3025) PEN0328 
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48 PRINT 7010 PEN0505 

7010 FORMAT ( 1H1////78H AIRCRAFT BEGINNING ELEVATION IS TO LOW TO AVOID PEN0506 

^TERRAIN. THIS RUN TERMINATED. ) PEN0507 
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Appendix III 



Inputs and Outputs of PENAIR 

Data cards for the PENAIR program should be prepared in the following manner 



Data Number 

Group of Cards 

1 1 



Variable Columns 

Name of Card 

NRUN 1-10 



Description 
of Variable 

A right justified integer, 
the number of runs to be 
performed . 



2 



KMAX 1-10 



XMIN 11-20 



YMIN 21-30 



XMAX 31-40 



A right justified integer, 
the number of squares used 
to approximate the terrain, 
KMAX s; 12 o 

A right justified floating 
point number, the x co- 
ordinate of the left edge of 
the terrain . All coordinates 
are in meters , that is the 
distance between the two 
x coordinates 72455 and 
72454 is one meter. 

A right justified floating 
point number, the y co- 
ordinate of the bottom edge 
of the terrain . 

The x coordinate of the right 
edge of the terrain. 
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I 






Date Number Variable Columns 

Group of Cards Name of Card 

YMAX 41-50 



TI 51-60 



TPD 61-70 



Description 
of Variable 

The y coordinate of the 
top edge of the terrain. 

A right justified floating 
point number, the time 
interval of the game 
(seconds) , TI > 0 . 

A right justified floating 
point number, the dimen- 
sion of the map square 
approximated by one poly- 
nomial (meters) . 



3 



4 



1 X 1-10 

Y 11-20 

Z 21-30 

X2 31-40 

Y2 41-50 



1 GMAX 1-10 



RC 11-20 



Right justified floating 
point numbers, the begin- 
ning coordinates of the 
aircraft. If IMARK = 3„4, or 
6,. y must equal YMIN. 

The x and y coordinates of 
the desired destination of 
the aircraft . If IMARK = 7 
or 8, this card may be blank. 

A right justified floating 
point number, the maximum 
positive G force to be used 
by the aircraft. 

A right justified floating 
point number, the aircraft’s 
maximum rate of climb 
(feet/minute) 
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. 



\ 










Data Number Variable Columns 

Group of Cards Name of Card 



Description 
of Variable 



RD 21-30 A right justified floating 

point number, the aircraft's 
maximum rate of descent 
(feet/minute) . 

GS 31-40 A right justified floating 

point number, the original 
airspeed of the aircraft. 



5 



IMAX 1-10 A right justified integer, 

the number of types of 
weapons, IMAX £ 10 . . 



6 



WPNAWPNB 1-16 
JMAX 21-25 



TA 26-30 



FT 31-35 



The name of the weapon, 
e . g . , 50 Cal . MG 
A right justified integer, 
the number of weapons of 
this type, JMAX ^ 20 . 

A right justified integer, 
the time required to acquire 
the target by this weapon 
(seconds) , TA s 0 . 

A right justified integer, 
the firing time of this 
weapon before reloading 
is required (seconds) , 

FT > 0. 
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. 









Data 


Number 


Variable 


Columns 


Description 


Group 


of Cards 


Name 


of Card 


of Variable 






RLT 


36-40 


A right justified integer. 










the time to reload this 










weapon (seconds) RLT £ 0 . 






RMAX 


41-50 


A right justified floating 
point number, the maximum 
range of this weapon(meters) 
RMAX > 0 . 






RMIN 


51-60 


A right justified floating 
point number, the minimum 
range of this weapon(meters) . 
RMIN ;> 0 . 


7 


2 


P 


1-8 


The probability of kill by 








9-16 


this weapon in one time 










interval at ranges of 








72-80 


100, 200, 300, . . . , 1000 










meters . 








1-8 


The second card of this 








9-16 


group contains the kill 










probabilities at ranges of 








72-80 


1100 to 2000 meters in 100 










meter intervals . 


8 


JMAX 


NW 


1-2 


The number assigned to the 










weapon, 1 £ NW £ 99 . 






XW 


11-20 


The x, y and z coordinates 






YW 


21-30 


of this weapon, ZW may be 






ZW 


31-40 


blank since PENAIR will com- 










pute the weapon's elevation. 
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NOTE: Data groups 6 , 7 and 8 should be repeated IMAX times . 



Data 


Number Variable 


Columns 


Description 


Group 


of Cards 


Name 


of Card 


of Variable 


9 


1 


MMAX 


1-10 


A right justified integer. 










the degree of x in the poly- 










nomial that approximates 










the first map square. 






NMAX 


11-20 


A right justified integer. 










the degree of y in the 










polynomial that approximates 










the first map square 


10 


As many cards 


PHI 


1-80 


The coefficients of the 




as necessary 






first polynomial appear on 



this and the following cards 
of this data group in a 
4E20.10 field as obtained 
from the terrain simulation 
program . 



NOTE: Data groups 9 and 10 should be repeated KMAX times . 



11 4 AFIELD 1-16 These four cards contain 

the four lines of the title 
to be printed at the top of 
the output. 



80 
















i 






' 

i 

: 



Data 


Number 


Variable 


Columns 


Description 


Group 


of Cards 


Name 


of Cards 


of Variable 


12 


1 


IMARK 


10 


The*value of this flag, deter 



mines the flightpath to be flown. 
IMARK = 1, upon receiving fire, 
fly straight and level , 

= 2, climb straight ahead, 

= 3, turn away, level, 

= 4, turn away and climb, 

= 5, dive and fly nap of 
the earth, 

= 6, turn away and fly 
nap of the earth, 

= 7, the flight path is 
completely pre- 
determined and read 
in as data group 13 , 

= 8, the x and y coordi- 
nates are read in 
as data group 13 , 
and a nap of the earth 
flightpath is computed,, 

IPUNCH 20 IPUNCH = 1, the flightpath will 

will be punched on 
cards . 

= 0, the flightpath will 
not be punched. 
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Data 


Number 


Variable 


Column 


Description 


Group 


of Cards 


Name 


of Card 


of Variable 






GSNEW 


21-30 


A right justified floating 



point number, the ground 
speed to be attained by the 
aircraft after receiving fire . 



NOTE: Data groups 11 and 12 will be repeated NRUN times. 



As many 
cards as 
necessary 
but not 


NPATH 


6-7 


Any positive integer, used 
to identify the flightpath. 


greater 


XT 


18-26 


Right justified floating 


than 300 . 


YT 


27-35 


point numbers , each of the 




ZH 


36-44 


cards of this group contain 
x, y, and z coordinates of 
the aircraft inconsecutive 
time intervals. 



NOTE: A blank card should be inserted at the end of this data group. 
This data group will be used only if IMARK = 7 or 8. If 
IMARK = 8, ZH may be blank. Data group 13 should follow 
its corresponding data groups 11 and 12 . 
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Output of the PENAIR Program 



A sample output from the PENAIR program is illustrated in Figures 
14 and 15. The input parameters of the weapons are printed with the 
weapon positions and elevations . The second page of the output (Figure 
15) is a summary of time in view and aircraft survival probability. The 
time in view listed here is in number of time intervals in view. If the 
time interval used in this example had been 0.5 second, the aircraft 
was in view of weapon number 49, located at (70080, 76940, 339), for 
28 seconds (56x0.5). The survival probability of the aircraft from this 
weapon was .9991. 
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SUMMARY C)F W LAP ON DATA 
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FIRST PAGE OF OUTPUT FROM THE PENAIR PROGRAM 
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FIGURE 15 

SECOND PAGE OF OUTPUT FROM THE PENAIR PROGRAM 



APPENDIX IV 



FORTRAN Listing of the TERRAIN Program 

TERRAIN consists of a main program and three overlays . The control 
cards for the COOP-Monitor of the CDC 1604 are included in this listing 
and are explained in Appendix V. The variable names used in the TERRAIN 
program are defined in the comments appearing in the FORTRAN listing. 

The variable names of input and output parameters are explained in Appendix 
V„ 
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DENC=DENC+ Q(M-1,J)**2 TER0282 
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FLMD=MD TER0495 
FLND=ND TER0496 
IF (IB) 51,50,51 TER0497 
10 = -1 TER0498 
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END TER0711 
SUBROUTINE TRACE ( AM »MY »NY ) TER0712 
DIMENSION AM ( MY ♦ NY ) > REC ( 900 ) > X(1800), Y(1800) TER0713 
DIMENSION I PT ( 3 » 3 ) » I NX ( 8 ) , I NY ( 8 ) TER0714 
DIMENSION DUMDUM ( 3 84 ) TER0715 
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I S= I S-8 TER0747 

I DX = I NX ( IS) ‘ TER07<+8 

I DY = I NY (IS) TER07A9 

IX2=IX+IDX TER0750 

IY2=IY+IDY TER0751 
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CALL CALC ( AM * MY * NY ) TER0783 

I X = I X+ I DX TER0784 

GO TO 6 TER0785 

I Y= I Y+ I DY TER0786 
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Z= I Y TER0819 

IY2=IY+IDY TER0820 

DY=IDY TER0821 

41 Y ( N ) = ( (AM( I X » IY)-CV)/(AM( IX»IY)-AM(IX* I Y 2 ) ) ) *DY + Z TER0822 

RETURN TER0823 
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241 SAVE1 = X ( NUMPTS + 1) TER0999 

SAVE2 = Y ( NUMPTS + 1) TER1000 

X ( NUMPTS + 1) = X(NUMPTS) TERlOOl 

Y ( NUMP T S + 1) = Y(NUMPTS) TER1002 

242 INUM = (NUMPTS + 1 ) / 2 TER1003 
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ITOP = JH - IY100 TER1035 
I RHT = LH - 1X100 TER1036 
NEXT 1 = ITOP T ER 10 37 
NEXT2 a I RUT TER1038 
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JGR I D ( J+4 ) = ICURV5 TER1071 

CALL IWRITE ( JGR I D » IDUMMY, JGRID(J + A)) TER1072 

I F ( IDUMMY) 5000*9025, 5000 TER1073 

SET UP RETURN. TER1074 

9025 IF (MODCURV-3) 87,278,87 TER1075 
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DIMENSION A ( 1800) TER1107 
F = A ( 1 ) TER1108 
DO 1 11=2, KK TER1109 
AA=A (II) TER1110 
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APPENDIX V 



Inputs and Outputs of the Terrain Simulation Program 

The following information is presented for potential users of this pro- 
gram: 

A. Tape units 

(1) The COOP control cards are shown in the FORTRAN listing (See 
Appendix IV.) A dash (-) in column one indicates a 7,9 punch, while an 
open parenthesis (() indicates an 11,0,7,9 punch,, 

(2) Logical tape unit 49 is the graph tape and 45 refers to the type- 
writer. 

(3) Scratch tapes 56 and 57 are necessary and therefore a total of 
eight tape units are required 

B. Input Formats 

The preparation of the input data cards is explained on the following pages 
of this appendix. 
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Data Number Variable Columns Description 

-Group of Cards Name of Card of Variable 



1 2 ITITLE 1-48 These two cards contain the 

two line title desired for 
the first page of the output 
and the bottom of the plotted 
map. 



2 



3 



1 NL 1-5 A right justified integer, 

the number of contour lines 
requiring a label. 



NL II 



LABL(Il) 



1-5 A right justified integer, 

the number of the contour 
level to be labeled , 

1 s II <: 50 

7-10 A left justified integer, 

the elevation of contour 
level II. 



4 1 NMAX 1-5 

MMAX 6-10 

NEWCTL 11-15 

is desired. 

= 1 , an output tape 

will be generated 
for the CalComp 
165 plotter. 



A right justified integer, 
the maximum degree of x to 
be considered 
A right justified integer, 
the maximum degree of y to 
be considered . 

A right justified integer, 
NEWCTL = 0 , no contour map 
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Data 

Group 

5 



6 



Number Variable Columns Description 

of Cards Name of Cards of Variable 



1 S 1-10 



A 11-20 



A right justified floating 
point number, the spacing 
desired between grid lines . 

S must be chosen so that 
neither the number of vertical 
nor horizontal grid lines 
exceed 45. 

A right justified floating 
point number, the linear 
dimension of the map square. 



1 



XTMIN 1-10 



XTMAX 11-20 



CND 21-30 



A right justified floating 
point number, the elevation 
of the lowest contour in 
the map square. 

A right justified floating 
point number, the maximum 
elevation in the map square. 
A right justified floating 
point number, the contour 
interval . The first contour 
is the value C = XTMIN. 



1 GX 1-2 The first two digits of a 

five digit floating point 
number, GX * L000 = the x 
coordinate of the lower 
left comer of the map square. 
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Data 

Group 



Number 


Variable 


Columns 


Description 


of Cards 


Name 


of Cards 


of Variation 




GY 


3-4 


The first two digits of a five 
digit floating point number , 
GY * 1000 = the y coordinate 
of the lower left comer of the 
map square. 




XP(1) 


11-13 


A right justified floating point 
number, the x coordinate of 
the lower left comer of the map 
square measured on an arbitrary 
scale. XP (1) usually should 
equal 0000 . 




YP(1) 


16-20 


The y coordinate of the lower left 
comer measured on an arbitrary 
scale. YP(1) usually should equal 
0000. 




XP(2) 


21-25 


The x coordinate of the lower 
right corner. 




YP(2) 


26-30 


The y coordinate of the lower 
right comer. 




XP (3) 


31-35 


The x coordinate of the upper 
right comer. 




XP(3) 


36-40 


The y coordinate of the upper 
right comer. 




XP(4) 


41-45 


The x coordinate of the upper 
left comer. 




XP(4) 


46-50 


The y coordinate of the upper 
left comer. 
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Data 

Group 

8 



NOTE: 



Number Variable Columns Description 

of Cards Name of Card of Variable 



As many 
cards as 




1-4 


necessary 


FZ 


5-8 



NC 9-1)0 



X(l) 11-15 



Y(l) 16-20 



21-80 



Any identifying code desired 
Col. 1-4 are not read. 

A right justified floating 
point number, FZ * 10 = 
the elevation (feet) whose 
coordinates appear on this 
card . 

A sequence number starting 
at 01 and increasing by one 
on each card following whose 
coordinates are those of the 
same contour line . 

A right justified floating 
point number, the x coordi- 
nate of a point on the contour 
line. 

A right justified floating 
point number, the coordinate 
of a point on the contour line. 

F 5.0 fields of additional 
points on the contour line 
in the same manner as 
Col. 11-20. 



One blank card is inserted after this group of data cards (eighth 
data group) for each contour line. One additional blank card is 
inserted if a new seventh data group is read. 
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Data 


Number 


Variable 


Columns 


Description 


Group 


of Cards 


Name 


of Card 


of Variable 


9 


1 




7-9 


The figure 999 should appear 
in these colums to signify 
that this is the last card 
of the data deck . 



NOTE: There should be no blank cards between the last eighth data group 
and the ninth data group. 
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Example: 

It is desired that grid square 7175 of the U. S Army Hunter-Liggett 
Military Reservation be approximated by a polynomial. Figure 17 illus- 
trates the data cards necessary 
The words 

HLMR CLOSE TERRAIN 
GRID SQUARE 7175 

will appear at the top of the first page of the printed output and at the 
bottom of the plotted map The minimum and maximum elevation appear- 
ing in this square are 900 and 1200 feet respectively. The contour interval 
is 20 feet. The maximum degree of X and Y to be considered is 20. The 
grid interval will be 50 meters and the square is 1000 meters on a side. 

Figure 16 illustrates how the original map was enlarged to read the 
contour coordinates . 



(001 2,3485) 



( 0000 , 0000 ) 



(3455, 3440) 



(3452, 0008) 



FIGURE 16 



DISTORTION CAUSED BY PROJECTION 
OF A MAP GRID SQUARE 



Distortion caused by projection is apparent in Figure 16 and is cor- 
rected by the first overlay. The corner coordinates are arbitrary and all 
contour coordinates are relative to them. If the points are measured by 
hand, the actual corner coordinates may be used and in this case would 
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HLMR CLOSE TERRAIN 
GRID SQUARE 7175 
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FIGURE 17 



A SAMPLE DATA DECK FOR THE TERRAIN PROGRAM 
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be (71000.. 7 5000) , (72000 , 75000) , (72000, 76000), (71000, 76000). The 
contour coordinates would then be the actual terrain coordinates . 

In this example the grid square was reprojected a second time with new 
corner coordinates of (0000, 0000), (3444,-0020), (3441, 3467), (-0004,3475). 

Terrain has three types of output, a) printed, b) punched cards, and 
c) a plotted map (See Figure 18) . The printed output consists of 1) the x and 
y coordinates of the grid to be used, 2) the input elevations at points on this 
grid as interpolated by Overlay 1,3) the elevations of the same points as 
computed using the approximating polynomial, 4) a table of differences be- 
tween the actual and computed elevations of the grid, 5) the computed ele- 
vations to be used in Overlay 3, 6) the coefficients of the approximating 
polynomial, and 7) the actual mean elevation squared, computed mean 
elevation squared and the maximum differences between the actual and com- 
puted elevations . Although accuracy is not determined in general , if the 
actual mean elevation squared and computed mean elevation squared are 
equal, then computation is exact. In practice they will not be equal due 
to the imprecise nature of calculation. A wide discrepancy indicates ex- 
cessive errors in calculation. 

The first card of the punched card output contains the degree of x and 
y in a 2110 format. The remaining cards contain the coefficients of the 
approximating polynomial in a 4E20.10 format. 
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tail 033 • 026 01? 003 001 



eai 






817 






ejj 



am 




HLMR CLOSE TERRAIN 
GRID SQUARE 7174 



FIGURE 18 

A SAMPLE OF A PLOTTED MAP 
FOR THE TERRAIN PROGRAM 
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